Splitting and merging an elongated Bose-Einstein condensate at finite temperature 
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We analyze coherence effects during the sphtting of a quasi one- dimensional condensate into 
two spatially separated ones and their subsequent merging into a single condensate. Our analysis 
takes into account finite-temperature effects, where phase fluctuations play an important role. We 
show that, at zero- temperature, the two split condensates can be merged into a single one with a 
negligible phase difference. By increasing temperature to a flnite value below the critical point for 
condensation (Tc), i.e., < T /Tc < 1, a considerable enhancement of phase and density fluctuations 
appears during the process of splitting and merging. Our results show that if the process of splitting 
and merging is sufficiently adiabatic, the whole process is quite insensitive to phase fiuctuations and 
even at high temperatures, a single condensate can be produced. 

PACS numbers: 03.75.Lm,03.75.Nt,05.30. Jp 
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INTRODUCTION 

The experimental exploration of Bose-Einstein conden- 
sates (BECs) in many different systems such as alkali 
metal gases 0,|3iM, hydrogen meta-stable helium 
and molecules 0, Xh, i^] has led to a surge of interest in 
manipulating ultracold atomic samples under very dif- 
ferent circumstances. One of the initial motivations for 
such development was and is the prospect of creating a 
continuous and coherent atomic beam 0, liol IllJ. which 
is the atomic analogy of the optical laser beam [l^l . 

Among the major challenges in achieving a continuous 
atom laser is how to overcome the difficulty due to the 
restrictive cooling conditions for continuously condens- 
ing the atomic gas. Spatial separation of the evaporative 
cooling from the destructive laser cooling may play a role 
in alleviating these challenges 0|. Recently, a continu- 
ous BEG source was achieved by periodically replenishing 
a reservoir with condensates There, optical tweezers 
were used to transport sodium condensates from where 
they were created to the reservoir. Such procedure allows 
one to continuously replenish the reservoir which contains 
on average more than 10^ atoms. Bringing a new con- 
densate to the reservoir each time implies the merging of 
condensates with different relative phases since each con- 
densates is created with a different phase. The merging of 
two condensates with different phases poses a limitation 
on the coherence of the process. 

Recently, interference effects in the merging of 30 
uncorrelated Bose-Einstein condensates released from 
a one-dimensional optical lattice have been discussed 
in whereas coherent splitting of BECs by deform- 
ing an optical single- well into a double-well potential for 
a trapped atom interferometer have been addressed in 
Refs. 0, 0|. Very recently, Schumm et al. [l^ has 
demonstrated a coherent beam splitter on an atom chip 



by splitting the condensate in a double- well potential and 
merging it again. They have demonstrated phase preser- 
vation in this process, even in the case when the split 
BECs are far enough apart to inhibit tunnel coupling. 

In this paper, we analyze the axial splitting of a very 
much elon gat ed cigar-shape condensate into two conden- 
sates 0, |3| and their subsequent merging along the 
axial direction. Our analysis includes finite-temperature 
effects. In other words, phase fluctuations arising due to 
temperature are taken into account during the whole pro- 
cess: splitting and merging. We observe that as long as 
the process of splitting and merging is adiabatic enough, 
both the split and the merged condensates, even at rela- 
tively large temperatures, do survive the effects of these 
fluctuations. 

Low-dimensional quantum gases exhibit very fascinat- 
ing properties and have attracted a lot of interest, both 
theoretically and experimentally [13 El El El IM HI 
E3- It is known that low-dimensional [one- (ID) and 
two- (2D) dimensional] quantum gases differ qualitatively 
from their 3D counterparts [UlilElEi- From a the- 
oretical point of view, the use of a mean-field theory to 
describe a low-dimensional quantum gas is severely re- 
stricted. A widely used criterion to apply a mean-field 
approach is to demand that the average distance between 
particles, d, is clearly smaller than the correlation length 
of the condensate Ic = h/ ^mng where m, g, and n de- 
note the mass, the interaction coupling, and the density, 
respectively. In three dimensions, the above condition 
leads to Ic/d cx n~^/^ and is well satisfied for small den- 
sities, and the description of the system at T = with a 
mean-field Gross-Pitaevskii equation is fully satisfactory. 
In the one-dimensional case, however, this ratio behaves 
as Ic/d oc v}/'^ and this fact changes drastically the range 
of applicability of a mean-field treatment. A careful anal- 
ysis of the validity of a mean-field treatment in the ID 
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case [30| leads to the following condition: 



where N is the number of condensed atoms, — 
y^h/{muJz) and a± — ^ hj {muj±) are the axial and ra- 
dial oscillator lengths, respectively, and ujz and uj± are 
the angular frequencies in the corresponding directions. 
Thus, in ID, contrary to the 3D case, when the density 
decreases the gas gradually becomes strongly correlated, 
acquires a fermionic character, and enters into the so- 
called Tonks-Girardeau regime j^IUH^Q- Experi- 
mental demonstration of a Tonks gas has been recently 
achieved 0. 

The possibility of generating low-dimensional bosonic 
gases raises the question of the effects of quantum fluctua- 
tions. In an untrapped ID Bose system these fluctuations 
destroy finite- as well as zero-temperature condensation. 
For trapped Bose gases, the situation is quite different: 
for noninteracting bosons in such a trap the finite spac- 
ing between the lowest and the next energy level allows 
for the occurrence of ID Bose-Einstein condensation even 
at finite-temperatures as stipulated in Refs. 113, In 
such a case the absence of gapless excitations indicates 
that the BEG will not be destroyed immediately as in- 
teractions between bosons are turned on. 

In ID trapping geometries, long-wavelength density 
and phase fluctuations lead to a new intermediate state 
between a condensate and a noncondensed system, which 
is commonly referred to as a quasicondensate. In quasi- 
condensates, the phase of the condensate is only coher- 
ent over a finite distance that is smaller than the sys- 
tem size. In other words, the phase coherence length is 
smaller than the axial size of the sample. To understand 
the nature of quasicondensates at finite-temperature, one 
has to analyze the behavior of the single particle corre- 
lation function by calculating the fluctuations of phase 
and density as has been done by Petrov et al. . There 
it is shown that for temperatures below the degeneracy 
temperature, the condensate's phase indeed fluctuates, 
but fluctuations of the density are still highly suppressed. 
This character of thermal fluctuations is also present in 
highly elongated 3D gases l3a , and has been recently 
observed experimentally in [MIsTlliiliiliollilt. 

Quasi ID geometries can be accessible in magnetic 
traps with a cylindrically symmetric harmonic potential 
along the axial direction that have transverse frequencies 
a;_L much larger than the axial one lo^. In such configura- 
tions the resulting condensate looks like a cigar aligned 
along the symmetry or z axis. With current technology, 
condensates with aspect ratio A — uo±^/ujz as large as 1000 
are achievable. 

In this paper, we study the process of splitting and 
subsequent merging of an elongated condensate confined 
in a ID geometry both at T = 0, i.e., when the conden- 
sate has a well-defined phase, and at finite-temperature 



(T > 0), in the quasicondensate regime. For finite T, 
we analyze the process of splitting and merging for a 
wide range of temperatures, i.e., < T < Td, where 
Tfji — Tdfiujz/ii being the chemical potential) corre- 
sponds to the characteristic temperature above which a 
true-condensate turns into a quasicondensate in which 
phase fluctuations begin to play a role. On the other 
hand Td = Nhuz/kB is the ID degeneracy tempera- 
ture |2l|. The transition, or crossover, between the dif- 
ferent regimes for the ID degenerate interacting bosonic 
gas, i.e., true-condensate, quasi-condensate, and Tonks 
gas is smooth. Thus, in the regime of quasi-condensation 
density fluctuations are relatively suppressed while phase 
fluctuations are enhanced. By keeping all parameters 
fixed and reducing the number of atoms, phase fluctu- 
ations become more and more pronounced, mean-field 
theory fails, and the gas enters into the strongly corre- 
lated regime or Tonks gas 

In our case, the splitting of the condensate is achieved 
by means of a double-well potential grown adiabatically 
on top of a trapping harmonic potential. By adiabati- 
cally switching off the double-well potential, we merge 
these condensates into a single one (Merged BEG). We 
would like to stress that we use the mean-field Gross- 
Pitaevskii equation (GPE) throughout this paper. It is 
worth pointing out that although the Gross-Pitaevskii 
equation [4^ describes properly coherent evolution of 
the atomic mean-field at T = 0, it can also be used 
to solve time evolution at finite-temperature in a rela- 
tively straightforward manner. It might look like that 
the mean-field method of the GPE allows to make state- 
ments only about first order coherence. But, as it is 
well known, a closely inspection reveals that the GPE 
contains as such classical Bogoliubov-de Gennes (BdG) 
equations, i.e. equations describing small fluctuations 
around a given solution of the GPE. Since we are here 
interested in the study of relatively high temperatures 
T > T^, classical description of fluctuations is fully ap- 
propriate. It turns out that one can simulate tempera- 
ture effects by adding fluctuations to the ground state 
solution of the GPE at T = in a way which mim- 
ics thermal fluctuations. At low temperatures T < Tc, 
where Tc = N / h\{2N)hwz/kB [ISl, this can be generally 
done by identifying phonon (quasi-particle) modes, i.e. 
eigensolutions of the BdG equations. The fluctuations 
are expressed thus as sums over the quasi-particles with 
amplitudes taken from Monte Garlo sampling and corre- 
sponding to the thermal (Boltzmann) populations of the 
quasi-particle modes. Such a method is used in this pa- 
per, with the additional simplification that for quasi-lD 
situations only the phase fluctuations are relevant. 

Note, that once we add the fluctuations initially at 
t = 0, and as long as they remain small in the course of 
evolution, they will propagate in time with a very good 
approximation as appropriate solutions of the time de- 
pendent BdG equations. Note also, that our approach 
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allows in principle to obtain information about all cor- 
relation functions: either in the form of an average over 
Monte Carlo realizations of the initial fluctuations, or in 
the form of time averages, due to the (expected) ergodic 
character of the evolution. In some cases, even averag- 
ing over the initial data is not necessary: the results for 
different realizations are so similar that looking for few 
single cases allows one to draw conclusions about "co- 
herence" . By coherence we mean here, that the splitting 
process of an initial ID BEC into two spatially sepa- 
rated ID BECs occurs with a well-defined relative phase 
between them. Their subsequent merging into a single 
merged BEC, when the process is fully coherent, should 
result in a single condensate with a well-defined relative 
phase with respect to the initial one. If this is the case, 
there is a perfect overlap between the densities of the 
initial and merged condensates. 

It is worth stressing that our approach is a simpli- 
fied version of the classical field methods used by sev- 
eral groups 

[UlillillillillillillillH^ 

In their 

approach some emphasis was put on the explanation of 
the temperature concept, and a certain model of finite- 
temperature effects has been studied. Particularly in- 
teresting is here the possibility of extracting higher order 
correlations from a single shot measurements |H3,IESIH3' 
which in our case corresponds to a single realizations of 
the initial fluctuations. 



DESCRIPTION OF THE MODEL 

Here we consider ^^Rb condensate with N = 1.2 x 10^ 
atoms confined in an harmonic trap with frequencies 
u}± = 2tt X 715 Hz and = 27r x 5 Hz. For such pa- 
rameters, the system is in the ID Thomas-Fermi regime 
{n ^ hojz) along the axial direction. The ID chemical 
potential is given by 



3 N m gii)\/h/mujz 
4\/2 



2/3 



(2) 



where q m = 2hujA_as is the effective ID coupling 
strength jSSj . For our parameters, transverse excitations 
are suppressed ~ fioJ±), and the dynamics of such a 
Bose gas can be described by the usual mean-field GPE 
in ID 
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+ VT{z,t) 



dt V 2 m dz^ 

-f 5iDiV|vI/(z,i)|2)M/(z,i), 



(3) 



where ^{z, t) is the mean-field order parameter, or in 
other words the condensate wave function. The potential 
term Vriz^t) includes both the magnetic trap and the 



double-well potential as described below: 

VT{z,t) = T4rap(z) + K>p(z,i) 

= \mujlz'^ + S{t)Vo cos(fc/ zf. 



(4) 
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FIG. 1: (Color online) Illustration of a double- well potential 
for splitting and merging of an elongated ID BEC: a strong 
optical potential and an axial trapping potential are combined 
for creating a double- well potential which is used for splitting 
a ID condensate, when switched on adiabatically, into two 
spatially separated symmetrical ID condensates at the center 
of each well, and then merge them into one when switched off 
adiabatically again. 



The potential used to split the condensate into two spa- 
tially separated condensates, I4p(z, i), is switched on and 
off adiabatically by means of a time-dependent function 
S{t). The maximum depth of this potential is Vb = 2.2 x 
IQ'^Er in terms of the recoil energy Er = fi^kf /2rn and 
ki — 27r/ A/ . To achieve spatial separation of the split con- 
densates one has to require that the distance between the 
two wells is, at least, of the order of the Thomas-Fermi ra- 
dius (i?TF-) This radius is given by R-jf = 2 ^ / {rriLol) , 
so that ki — 7r/i?TF- For N ~ 10^ atoms, the Thomas- 
Fermi radius is i?TF — 88/im. This is in agreement with 
the experimental results in Refs. 0, 0| . 

The time dependent function S{t) controls the switch- 
ing on and off of the double- well potential and hence the 
overall splitting and merging process of the condensates. 
We define this function as 



S{t) 



0, 

sin(f: 
1, 

cos 



for t < t„vo and t > t„ 



for tr. 



<t<t 



spli 



for tspi <t<tc 



(5) 



In this equation Q, ^cvo is the time duration required 
for evolving the ID BEC in real time before the splitting 
process begins. It extends for 10ms. Just immediately 
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after, the splitting process begins and continues for a 
time interval tgpi- At the end of tspi, the function S(t) 
attains a maximum value of unity and remains constant 
for a time interval icon- During this time interval, two 
spatially separated ID BECs are created. At the begin- 
ning of the time span i„icr, merging of the condensates 
starts by switching off the doublc-wcU potential. This 
process continues until the two ID BECs merge together. 
A complete merging is only possible when S{t) becomes 
finally zero, i.e., when the optical potential is completely 
switched off. In this case the atoms remain only under 
the influence of the trapping potential. 

To ensure coherence during the process of the switch- 
ing on and off of the optical potential, the raising of 
the double-well has to be slow enough to avoid excita- 
tions and to allow for quantum tunneling between the 
two wells. Notice that the relevant time scale of excita- 
tions is given by the inverse of the frequency of the trap, 

in our case tsys = Ztt/Wz = 200 mS. If fspl,imer ^ ^sys 

we expect coherent splitting and merging. On the other 
hand, for tsp\, imcr of the order of tgys or less, the process 
creates more and more excitations that cause incoher- 
ence. The faster the double-well potential switches on 
and off the stronger the excitations are. To check our 
claims we have carried out further numerical simulations 
for several values of tspi and tmcr around tsys- In general, 
we have observed that a coherent process, at T = 0, is 
achieved for switching times of 400 ms (or larger). 

The full time dependency of S{t) necessary for the co- 
herent splitting and merging of the condensates based on 
Eq. (O is illustrated in Fig. 13 




FIG. 2: (Color online) Time dependence of the function S{t). 
For the whole process, we fixed a time evolution of tcvo ~ 10 
ms, a splitting time of t^pi — 440 ms, a time icon = 430 ms in 
which S(t) remains constant (unity), and a time of merging 
imcr ~ 530 ms and an additional time of 10 ms for allowing 
the merged condensate to make its final evolution into a single 
condensate. 



SPLITTING AND MERGING AT ZERO 
TEMPERATURE 

We calculate first the ground state of the system by 
evolving the GPE (|3J| with S{t) = in imaginary time. 



After the ground state has been found, we numerically 
solve Eq. (PJ with the pulse profile given by the time 
dependent-function S{t). 

The results of our simulations at temperature T = 0, 
are summarized in Fig. |21 There we display the conden- 
sates' density at three different times. First at < = 0, 
the initial condensate (Initial ID BEC) has a character- 
istic Thomas-Fermi profile with i?TF = 88/im. Then at 
t = 800 ms (corresponding to the ramping up of the 
double-well potential in 400 ms and keeping it constant 
during additional 400 ms), two spatially well separated 
identical condensates (Split BECs) appeared, centered 
at z/i?TF = ±0.5, each of them with a number of atoms 
Ns = N/2. Finally, we display the density of the conden- 
sate a.t t — 1400 ms. The merged condensate (Merged 
ID BEC) has exactly the same profile as the initial one 
(Initial ID BEC) and thus they cannot be distinguished 
in the figure, asserting that the process of splitting and 
merging is fully coherent. 




FIG. 3: (Color online) Adiabatic splitting and merging of 
an elongated condensate at T = 0. After coherently split- 
ting an initial ID condensate (Initial BEC) into two spa- 
tially separated ones (Split BECs), by adiabatically switch- 
ing on the double-well potential given in Fig. centered at 
z/Rtf ~ —0.5 and z/Rtf = 0.5, the optical potential is again 
switched off adiabatically. This leads to a coherent merging of 
the Split condensates into a single one (Merged ID BEC). The 
merged ID BEC overlaps on top of the initial ID BEC which 
is an indication of a completely coherent merging process. 

Coherence may be a prerequisite for further appli- 
cations such as in atom interferometry and quantum- 
information processing |54l| . 

Before we proceed to the case of finite-temperature, it 
may be relevant to address the case of a nonadiabatic 
splitting and merging process, i.e., when the double- well 
potential switches on and off too fast. We have carried 
out simulations for switching times as short as 20 ms. In 
such cases the splitting of the condensate (even at T = 0) 
becomes completely incoherent, there is no trace of phase 
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preservation, and the condensate is destroyed, as is shown 
in Fig. 2] For shorter sphtting times, tgpi < 20 ms, it is 





Initial BEC 

Incoherent Split 
^yicoherent Merging - 


/ 




/ J 


illlJiil| \ 












/ ' lili 


i;?' ■ i' 'i ' ' 1 11* 





FIG. 4: (Color online) Nonadiabatic splitting and merging 
of an elongated condensate at T = 0. The Initial BEC 
is split incoherently by considering a very small splitting 
time (tspi) into two picked spatially separated ones (Incoher- 
ent Split) by nonadiabatically switching on the double-well 
potential (Fig. 0. The optical potential is again switched 
off nonadiabatically which leads to a completely incoherent 
pattern (Incoherent Merging) due to larger excitations. 

not even possible to split the Initial BEC into two well 
spatially separated condensates. When it comes to the 
time of merging (imcr), a similar behavior is observed if 

^mcr — ^sys ■ 



SPLITTING AND MERGING AT FINITE 
TEMPERATURE 

So far we have considered only the case of T = where 
the initial condensate has a well-defined phase. In this 
section, following exactly the same approach as in the 
previous section, we investigate the effects that finite- 
temperature might have on the overall coherence during 
the splitting and merging process. 

Fluctuations of phase and density of a BEC are the 
result of thermal excitations, and appear usually at finite- 
temperature [s^ • In such a case the system includes 
some noncondensed or thermal particles and hence the 
total Bose field operator can be expressed as 



#(z,i) = 4'(z,t) -|-(5*(z,t). 



(6) 



where S^{z,t) describes the thermal depletion part. For 
a BEC in 3D trapping geometries, fluctuations of den- 
sity and phase are only important in a narrow tem- 
perature range near the BEC transition temperature 
Tc 113 ■ For ID systems, however, phase fluctuations are 
present at temperatures far below the degeneracy tem- 
perature. Phase fluctuations can be studied by solving 



the Bogoliubov-de Gennes equations describing elemen- 
tary excitations. Writing the quantum field operator as 
^(2) = y/niuiz)exp[i(t>{z)] where nijj{z) = \^{z,t)\'^ 
denotes the condensate density at T = [niY){z = 0) = 
/i/giD] the phase and density operators take, respectively, 
the following forms [sif: 



_^ [/+(.)«, +/-(z) a]] (7) 



and 



nioiz) = \/niB{z)^i [j^ % - /+ a] j , (8) 

where aj(aj) is the annihilation (creation) operator of 
the excitations with quantum number j and energy ej — 
huJz\/ j{j + l)/2, and = Uj ± vj, where Uj and vj 
denote the excitation functions determined by the BdG 
equations. More explicitly, the functions in a ID 
Thomas-Fermi regime and at finite-temperature take the 
form: 



(i + 1/2) /2/i . 

xPj{z/Rtf) 



{z/Rtv? 



±1/2 



(9) 



where Pj^z/Rtf) are Legendre polynomials. The phase 
coherence length, in terms of the Thomas-Fermi radius 
i?TF7 is expressed as = RtpTci ^/Thuz and character- 
izes the maximal distance between two phase-correlated 
points in the condensate. Phase fluctuations increase for 
large trap aspect ratios and small N 36] . 

Temperature is included at the level of the CPE 
by calculating flrst the density at T = in the presence 
of the magnetic trap only, and then mimicking flnite- 
temperature effects via the phase operator of Eq. (Cjl. In 
other words a phase is imprinted on the condensate wave 
function at this stage. To this aim, we calculate the Bose 
occupation Nj = [^e.'^i/^^B t) _ modes in the Bogoli- 
ubov approximation for fixed temperatures replacing the 
operators aj and by random complex variables aj and 
a*, respectively such that (|q!jP) = Nj [3fij |. 

Although the CPE © in this limit remains valid, the 
BdG equations become modified by the i nteg ration over 
the transverse profile of the condensate [33 ■ In effect 
the mode functions are given by Jacobi polynomials, 
whereas ej are given by a slightly different expression 
than in the pure ID case. We stress, however, that for 
the regimes of temperatures we consider, there will be 
no qualitative and practically no quantitative difference 
between the pure- and quasi- ID results. 

The dependence of the magnitude of the phase fluc- 
tuations of the condensate on the temperature is shown 
in Fig. 13 Even though we present here results calcu- 
lated for a flxed number of condensate atoms (-/V), it is 
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numerically verifiable that the magnitude of the phase 
fluctuations is inversely proportional to the square root 
of the number of atoms and hence to the peak density 
of the condensate as described by Eq. (Q. As can be 
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FIG. 5: (Color online) Enhancement of phase fluctuations 
with the increase in temperature, (a) At the top, we have very 
weak phase fluctuations at very low temperature (T/Tc — 
0.01). (b) With the rise in temperature {T/Tc = 0.1), the 
phase fluctuations begin to be enhanced, (c) Still at an inter- 
mediate but relatively high temperature (T/Tc ~ 0.4), the 
magnitude of phase fluctuations increases, (d) Finally at 
T/Tc ~ 0.8, which is near the critical point, stronger phase 
fluctuations are displayed. 

seen by inspecting the different plots (a)-(d) in Fig. 
the phase fluctuations get more and more enhanced with 
the increase of temperature. 

Having seen the enhancement of phase fluctuations 
with temperature, we proceed now to analyze the co- 
herence properties of splitting and merging of a ID BEC 
in the presence of these phase fluctuations for temper- 
atures in the range of < T < Tc. Our results for 
finite-temperature are summarized in Fig. |S1 

Using exactly the same approach that we employed for 
the adiabatic case at T = but now including temper- 
ature we study again the splitting and merging process 
for the same parameters. In this case, the density fluc- 
tuations, which are highly suppressed in elongated 3D 
condensates, are very pronounced in the ID density after 
time evolution, when the temperature increases from a 
very small value {T/Tc = 0.01) to a value near the crit- 
ical point for condensation {T/Tc — 0.8). At high tem- 
peratures, the density fluctuations get more and more en- 
hanced. However, notice that the effect of fluctuations is 
very similar on the single and on the merged condensates. 
On the other hand, fluctuations on the split condensates 
remain relatively small due to fact that the densities in 
the double wells are higher as shown in all plots of Fig. El 
This is in confirmation of the prediction of Eq. Q . The 
presence of fluctuations on the density profile is a conse- 
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FIG. 6: (Color online) Splitting and merging of ID BECs at 
finite-temperature, i.e., in the presence of phase and density 
fiuctuations. The four plots are the fluctuating densities at 
T/Tc = (I) 0.01, (II) 0.1, (III) 0.4 and (IV) 0.8. Each one 
in turn consists three plots with in it: curve a, the initial ID 
BEC (red solid line); curve b the split ID BECs (green dashed 
line); and curve c, the merged ID BEC (blue dotted line). 

quence of the BdG equations present in the GPE equa- 
tion. In spite of such fluctuations, the split and merged 
condensates present almost the same Thomas- Fermi den- 
sity profiles for any temperature T < Tc. 

From these observations, we conclude that in spite 
of the initial phase fluctuations in the quasicondensate 
regime at finite-tempcraturc, there is a preservation of 
"phase coherence length" during the splitting and merg- 
ing process, if adiabaticity is satisfied. 



SUMMARY 

Summarizing, we have discussed coherence effects in 
the splitting and merging of a ID ^^Rb BEC. This is done 
by creating two spatially separated condensates from an 
initial ID condensate by deforming the trapping poten- 
tial into a double-well potential. We have analyzed the 
case of zero temperature as well as finite-temperature in 
the so-called quasi-condensate regime, where the phase 
coherence length is smaller than the size of the system. 
At zero temperature and for a process adiabatic enough 
where the splitting and merging times are much larger 
than the characteristic time of the system given by the 
inverse of the trap frequency igys = Stt/wx, a coherent 
splitting followed by a coherent merging is achieved and 
there is a constant relative phase between the initial and 
the final merged condensates. On the contrary if the 
splitting and merging times are not larger than the rele- 
vant time scale, the split condensates acquire a random 
relative phases and merging is no longer possible. In 
this case, the system acquires large density and phase 
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fluctuations on a length scale comparable with the co- 
herence length Ic, and a description based on the GPE 
becomes invalid. In the case of finite-temperature our 
results show that even in the presence of strong phase 
fluctuations, if the process of splitting and merging is 
carried out fully adiabatically, the condensate preserves 
the Thomas-Fermi density profile and there is phase co- 
herence length preservation. In such cases, the merged 
condensate is a quasi-condensate with the same initial 
density profile as the initial condensate and with the same 
phase coherence length, although the "overall" phase of 
the quasi-condensate is not preserved. This situation oc- 
curs even at temperatures very near the critical value for 
condensation as long as the trapping potential remains in 
place. Our results may have a useful implication for ma- 
nipulating ID BECs at zero as well as finite-temperatures 
such as in atom lasers, interferometry, and solitons. In 
particular, our results agree qualitatively well with the 
recent measurements in Ref. llTl. 
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